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Abstract 

Viscosity plays a critical role in determining the stability of rotating neutron 
stars. We report the results of a calculation of the shear viscosity of /3 - stable 
matter, carried out using an effective interaction based on a state-of-the-art 
nucleon-nucleon potential and the formalism of correlated basis functions. 
Within our approach the equation of state, determining the proton fraction, 
and the nucleon-nucleon scattering probability are consistently obtained from 
the same dynamical model. The results show that, while the neutron contri- 
bution to the viscosity is always dominant, above nuclear saturation density 
the electron contribution becomes appreciable. 
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1. Introduction 

The quantitative description of transport properties of nuclear matter is 
relevant to the understanding of a variety of neutron star properties. Ther- 
mal conductivity is one of the driving factors of the cooling process, while 
electrical conductivity is relevant to the ohmic dissipation of magnetic fields 
in the star interior. In rotating stars, a crucial role is also played by vis- 
cosity, that determines the possible onset of the gravitational-wave driven 
instabilities first predicted by Chandrasekhar in the 1970s [H, 0|. 

Gravitational radiation (GR) is emitted when a non-radial oscillation 
mode of the star is excited by an internal or external perturbation. In non 
rotating neutron stars the emission of GR is a dissipative process, leading 
to the damping of the oscillation, while in rotating stars the effect can be 
quite different. At the end of the 1970s, Friedman and Schutz [3J proved 
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that, due to the mechanism discussed in Refs.[H, 0|, allperfect fluid rotating 
stars are in fact unstable. More recently, Andersson [4J and Friedman and 
Morsink [H| demonstrated that all the so called r modes, i.e. the oscillations of 
rotating stars whose restoring force is the Coriolis force, are driven unstable 
by GR in all perfect fluid stars. On the other hand, if matter in the star 
interior does not behave as a perfect fluid, dissipative processes, such as 
viscosity, can damp the modes responsible of the instability, or even suppress 
them completely. Hence, knowledge of the viscosity of neutron star matter 
is required to determine whether a mode is stable or unstable 

The main limitation of most analyses of the damping of neutron-star os- 
cillations [sl lies in the lack of consistency between the dynamical models 
used to obtain the equation of state (EOS), describing the equilibrium prop- 
erties of the star, and those employed to describe transport properties. In 
their seminal paper, Cutler and Lindblom [f| combined a variety of EOS, re- 
sulting from different theoretical approaches, with the pioneering estimates 
of the shear viscosity coefficient of neutron star matter obtained in the 1970s 
by Flowers and Itoh, who used the the Landau- Abrikosov-Khalatnikov 0, 0] 
formalism and the neutron- neutron collision probability estimated from the 
measured scattering phase shifts 0, • 

Nuclear many body theory provides a consistent framework to obtain the 
in medium nucleon-nucleon (NN) cross section and the transport coefficients 



of nuclear matter from realistic NN potentials, using either the G-matrix [11 



or the CBF 12J formalism. In both approaches one can define a well behaved 
effective interaction, suitable for use in standard perturbation theory in the 
Fermi gas basis and allowing; for a unified treatment of equilibrium and non 



equilibrium properties [ll|, [12|, [13|, [14], [15 



The CBF effective interaction has been employed to carry out a calcu 



lation of the shear viscosity of pure neutron matter [121 ] . However, a more 
realistic model of neutron star matter must allow for the presence of protons 
and electrons. As matter density increases, the electron chemical potential 
may also exceed the muon rest mass, making the appearance of muons ener- 
getically favorable. 



In this Letter, we discuss the generalization of the approach of Ref.[l2 
to the case of /^-stable matter consisting of neutrons, protons and electrons. 
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2. Formalism 



The application of the Abrikosov-Khalatnikov |8j formalism to the cal- 
culation of the shear viscosity of matter consisting of neutrons, protons and 
electrons in /^-equilibrium was first developed by Flowers and Itoh 0, [lo| . 
Within their approach, the NN scattering rate due to strong interactions 
is modeled using the measured free space cross section, thus neglecting all 
modifications caused by the presence of the nuclear medium. 

In Refs.jil, [loj, the calculation of the transport coefficients of /3-stable 
matter is carried through a straightforward generalization of the case of pure 
neutron matter. In a multicomponent system the Boltzmann-Landau equa- 
tion takes the form 

dna^dna de pa dn a de pa _ \^ j /-^ 
dt dr dp dp dr ^ a/3 ' 

where n a = n pa (r,t) denotes the distribution of quasiparticles of type a 
(a = n,p,e), carrying momentum p and energy e pa . The form of the 
collision term in the right hand side of the above equation clearly shows that, 
in principle, all binary collisions, involving both like and unlike quasiparticles, 
must be taken into account. 

The shear viscosity, defined as the coefficient of the momentum flux tensor 



appearing in the left hand side of Eq. ([!]), can be written as [10 



V = Vn + Vp + Ve , (2) 

the contribution associated with quasiparticles of type a being given by 0, 

IS EE El 



1 2 

Va = -p a mlv 2 Fa T a — —r C(£ aa ) . (3) 

5 ir 2 (l-£ aa ) 

In the above equation, p a , m* and VFa denote the density, effective mass and 
Fermi velocity, respectively, while the quasiparticle lifetime r a is given by 

4?r 4 



where T is the temperature, and 
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In Eq.flS]), W a p denotes the probability of collisions between quasiparticles 
of type a and (3. In the low temperature limit, underlying the Landau- 
Abrikosov-Khalatnikov approach, scattering processes can only involve quasi- 
particles carrying momenta close to the Fermi momentum. As a consequence, 
at fixed baryon density, p = p p + p p , and proton fraction x = p v j p, W a p only 
depends on two angular variables, 9 and <p, and the averages in Eqs.(jH) and 
(ED are defined as 

(F) = /S F(M) - (6) 

The quantities L°p and L a p are also functions of the angles 9 and (p. Their 
explicit expressions, as well as that of the factor C(£ aa ), are given in Refs. 0, 

m. 



3. Results 

For any given baryon density, the calculation of the shear viscosity re- 
quires the knowledge of the proton fraction x, determined by the conditions 
of /3-equilibrium and charge neutrality 

ftn Pp Pe i Pp Pe ; (7) 

where p a denotes the chemical potential of quasiparticles of type a. 

In this work, the proton and neutron chemical potentials have been com- 
puted within the Hartree-Fock approximation 

Pa = e a (p Fa ) , (8) 

where the Fermi momentum is given by pFa = (37r 2 p Q ,) 1//3 , using the single 



particle spectrum obtained from the CBF effective interaction of Ref . [12 
the Hartree-Fock approximation 

e a (p) = x a + p^Xp J d 3 x [{v eS ) D - (v eS ) E £(pF a x) e ip ' x ] J . (9) 

In the above equation, m denotes the nucleon mass, £(x) = 3 (sin x—x cos x) / x 3 
and the spin averaged direct and exchange matrix elements of the effective 
interaction are given by 

(v c r)d = 9 XI ( a ^l u eff|a/?) , (v cS )e = -z X {<*P\ v &\P a ) ■ (10) 
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Figure 1: Solid line: single particle spectrum evaluated using Eq.([9|) and the effective 
interaction of Ref. [l2j . The diamonds show the results of Ref. 21 1 , obtained using the 
FHNC approach and a realistic nuclear hamiltonian. 



The effective interaction v c g is based on a truncated version of the NN po- 



tential referred to as Argonne V\$ [18j], providing an excellent fit of deuteron 



properties and the full Nijmegen phase shift database. It also includes the 
effects of interactions involving three- and many-nucleon forces, described 
according to the approach originally proposed in Ref. 19j. The EOS of sym- 



metric nuclear matter and pure neutron matter obtained using the dynamical 
model of Ref. 12] turn out to be in fairly good agreement with the results of 
the state-of-the-art calculations of Ref. 20j, carried out using the full Argonne 
v is potential. 

In Fig. [T]the single particle spectrum of symmetric nuclear matter at equi 
librium density obtained from Eq. (Q is compared to the results of Ref. 21 



carried out within the Fermi Hyper-Netted Chain (FHNC) approach using a 
realistic nuclear hamiltonian. 

Figure [2] shows the baryon density dependence of the proton fraction 
resulting from the numerical solution of Eqs.fJTJ), carried out assuming that 
electrons can be described using the single particle spectrum of the relativistic 
Fermi gas. The effective masses m* appearing in Eqs.([3]) and (j3J) can be 
readily obtained from the single particle energies through 



1 de(p) 
p dp 



(11) 
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Figure 2: Proton fraction of charge neutral /3-stable matter consisting of neutrons, protons 
and electrons, obtained from the effective interaction of Ref.[l2j|, as a function of baryon 
density. 



The scattering probabilities employed in our calculations take into ac- 
count both strong and electromagnetic interaction. Nuclear interactions 

W pn have been described within the 



contributing to W n 
dynamical model of Ref.|12|], while the calculation of W ep =W pe , W e n=Wn 



Wja,. and W np - 



and the electromagnetic part of W pp have been carried following Refs. [9|, 110 



Figure [3] shows the temperature-independent quantities r] a T 2 , as well as 
rjT 2 , plotted clS db function of baryon density. It appears that the proton 
viscosity r] p is always very small, due to their low density and mobility. On 
the other hand, in spite of the fact that p p =p e , the electron viscosity rj e turns 
out to be much higher, as electrons are ultra-relativistic. As expected, the 
dominant contribution IS Tj n ^ clS the neutron fraction is larger than 90% over 
the whole range of baryon density. For p < p , p = 0.16 fm ~ 3 being the 
equilibrium density of symmetric nuclear matter, the total viscosity can be 
identified with the neutron contribution. The electron contribution becomes 
barely visible only at larger density. 

The main difference between the case of a mult i- component fluid and 
that of pure neutron matter, discussed in Ref. 12|, lies in the larger number of 



channels available in collision processes, which implies a shorter quasiparticle 
lifetime, leading in turn to a lower viscosity. This feature is illustrated by 
the diamonds of Fig. [31 showing the energy dependence of the viscosity of 
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Figure 3: Baryon density dependence of the quantity r] a T 2 for protons (dotted line), 
electrons (dot-dashed line) and neutrons (dashed line) in charge neutral /3-stable matter. 
The solid line corresponds to the total r]T 2 (see Eq.Q). For comparison, the diamonds 
show rjT 2 of pure neutron matter, corresponding to x p = 0. 



pure neutron matter, corresponding to x p = 0. 

To gauge the dependence of our results on the composition of matter, we 
have repeated the calculations using a larger proton fraction, obtained by 
multiplying the results shown in Fig. [2] by a factor two. Note that changing 
the proton fraction amounts to using a different model of nuclear dynamics. 
Larger proton fractions correspond to a larger contribution of the symmetry 
term to nuclear matter energy. 

As shown in Fig. HI in this case the electron viscosity exceeds neutron 
viscosity at values of baryon density p ~ 0.2 fm -3 , just above nuclear matter 
saturation density, corresponding to an electron fraction of ~ 8%. 



4. Conclusions 

The results of our work show that inclusion of medium effects on NN 
scattering leads to a sizable suppression of the collision probability, producing 
in turn an enhancement of the shear viscosity of nuclear matter. 

The role played by screening of the bare NN interaction, mainly due 



to short range correlations, was already pointed out in Ref. [12j for the 
case of pure neutron matter. Our study of charge-neutral /3-stable matter 
consisting of neutrons, protons and electrons, carried out using the same 
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Figure 4: Same as in Fig. [3J but with the proton (and electron) fractions of Fig. [5] 
increased by a factor two. 



effective interaction, suggests that, while the proton contribution to viscosity 
can be safely neglected, the electron contribution is appreciable. In the case 
of dynamical models predicting a larger proton fraction, typically x p > 8% 
at p ~ po, it may in fact become dominant at densities exceeding nuclear 
matter saturation density. Overall, due to the availability of a larger number 
of reaction channels, the viscosity of /^-stable matter turns out to be lower 
than that of pure neutron matter. 

All the above considerations are based on the tenet that neutrons, protons 
and electrons behave as normal Fermi liquid. In the region of density and 
temperature in which neutrons become superfluid, their contribution to the 
viscosity vanishes and the electron contribution takes over. 

In their study of the effect of viscosity on neutron star oscillations, the 
authors of Ref. p took into account the onset of superfluidity writing the 
viscosity in the form 

ri a = [l-e(j>,T)]ri + 6(p,T)r) e , (12) 

where 

CO T > T 

e(p,T)= \ ^ , (is) 



T c = T c (p) being the critical temperature, taken from Ref.(22|. It has to be 
pointed out, however, that a fully consistent analysis of the stability of ro- 
tating neutron stars requires that the EOS, determining the proton fraction, 
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the viscosity coefficient and the critical temperature be all obtained from 
the same dynamical model. The effective interaction approach appears to be 
ideally suited to pursue this project. 

Numerical calculations of the 1 So superffuid gap carried out using the 
effective interaction of Ref. 12j yield a critical temperature T c ~ 2 x 10 10 K 
at density p ~ 0.04 fm -3 , typical of the the neutron star inner crust 23J, in 



fairly good agreement with the rsults of Ref.|22j. The extension of this study 
to the case of pairing in 3 P2 states is currently being carried out. 

As a final remark, it is worth mentioning that a more realistic model of 
neutron star matter should include muons, whose appearance is likely to be 
energetically favored at densities above nuclear saturation density. However, 
compared to electrons, muons have larger mass, and therefore lower mobility. 
As a consequence, taking into account their contribution to the viscosity is 
not expected to significantly affect the conclusions of our work. 
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